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ABSTRACT 

Context. The increasing number of observed supernova events allows for finding ever more frequently the progenitor star in archive 
images. In a few cases, the progenitor star is a yellow supergiant star. The estimated position in the Hertzsprung-Russell diagram of 
these stars is not compatible with the theoretical tracks of classical single-star models. 

Aims. According to several authors, the mass-loss rates during the red supergiant phase could be underestimated. We study the impact 
of an increase in these mass-loss rates on the position of 12 to 15 Mq stars at the end of their nuclear lives, in order to reconcile the 
theoretical tracks with the observed yellow supergiant progenitors. 

Methods. We have performed calculations of 12 to 15 Mq rotating stellar models using the Geneva stellar evolution code. To account 
for the uncertainties in the mass-loss rates during the RSG phase, we increased the mass-loss rate of the star (between 3 and 10 times 
the standard one) during that phase and compared the evolution of stars undergoing such high mass-loss rates with models computed 
with the standard mass-loss prescription. 

Results. We show that the final position of the models in the Hertzsprung-Russell diagram depends on the mass loss they undergo 
during the red supergiant phase. With an increased mass-loss rate, we find that some models end their nuclear life at positions that 
are compatible with the observed position of several supernova progenitors. We conclude that an increased mass-loss rate (whose 
physical mechanism still needs to be clarified) allows single-star models to simultaneously reproduce the estimated position in the 
HRD of the YSG SN progenitors, as well as the SN type. 
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During the past decade, the increase in the detection number of 
supernova (SN) events, combined with the growing num ber of 
identifications of th e progenitor star in archive images (iLi et alJ 
2007; Smartt 2009), become an important test for stellar evolu- 
tions models. SNe are classifi ed in several cate gories, based on 
their light curve and spectra (Filippenko 1997). The two main 
groups are the type II SNe (exhibiting hydrogen lines in their 
spectrum) and type I (without such lines). Most core col lapse 
SNe are type IIP SNe (ISmartt et al.ll2009USmith et al.ll201 ll) . re- 
quiring a progenitor with massive hydrogen envelope that pro- 
duces the typical "plateau" in the light curve. The plateau is 
due to the photosphere moving inwards through the mass of the 
ejecta while the ejecta expands, leading to a constant radius pho- 
tosphere and a plateau in the lightcurve. Fewer core-collapse 
SNe are found to be type IIL or type lib, whose progenitor is 
thought to have a much less massive hydrogen envelope. These 
two types of SNe require progenitors that encounter stronger 
tha n usual mass loss to partially remove their hydrogen-rich lay- 
ers dHeger et al.ll2003blEldridge & ToutfcOOH) . Finally the type I 
core collapse supernovae are separated into two subtypes, type lb 
(with helium features in the spectrum) and Ic (without helium). 

Among recent observations of supernova progenitors there 
are some cases that are challenging to explain in the frame 
of the c lassical single - star e volution. The progenitor of SN 
2008ax (ICrockett et al.1 [2008) lies in the Hertzsprung-Russell 
diagram (HRD) between the two main areas where the clas- 
sical non-rotating stellar model tracks end: the red supergiants 



(RSG) branch, on the one hand (with final log(r eff /K) ~ 3.55), 
and the Wolf-Rayet (WR) location, on the other (with final 
log(r e jf/K) ~ 4.5 - 5). However, the position of the progen- 
itor star can be explained by a singl e star leading to a WN 
star at the end of the evolution (ICrockett et al.ll2008h . and re- 
cent theoretical tracks with rotation and increased mass-loss rate 
allow reproduction o f the estimated position of its progenitor 
(Ekstrom et al.l201 lb . Severa l other SNe arise from stars that are 
not R SG: type lib SN 1993 J (lAldering et al.l [l994| : iMaund et alJ 
12004 . t ype IIP SN 2008cn (lElias-Rosa et aljborm type IIL SN 
2009kr dFraser et alj|2010t lElias-Rosa et aUBOlOl) and type lib 
SN 201 ldh dMaundetal.1120111) . all have a yellow supergiant 
(YSG) progenitors^, with surface effective temperature in the 
range 3.6 < log (T s g/K) < 3.8. For most of these progenitors, 
the observational data could be compatible with a binary sys- 
tem, which could explain the unusu al position of these stars in 
the HRD just before th eir explosion (Podsiadlowski et al. 1993; 
lEhas-Rosa et"ai]|2009l) . This is particularly true for SN 1993J, 
which is most probably a binary, and where strong observa- 
tional constraints exist about the secondary, and where the pro- 
genitor can be well explai ned by a binary stellar model (see 
Stancliffe & Eldridge 2009, and references therein). However, 
the recent obser vation of SN 201 ld h seems to indicate a single- 
star progenitor dMaund et al.ll2.01 lb . Single-star evolution mod- 
els should thus be able to explain at least this last peculiar case, 



1 The progenitor of SN 1993 is classified as a K supergiant 
fMaund et alj|2004l) . However, its log(T eff /K) = 3.63 makes it a "hot" 
RSG instead of a YSG. 
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Fig. 1. Left panel: HRD for our set of models. The initial mass and mass-loss prescription is indicated near the tracks. The circles 
indicate the end points of th e simulations . Right panel: Zoom on the region wh ere th e end points o f the t racks are situated. The 
observational data are from iMaund et all (120041) . lElias-Rosa et all (120091 feOlOl) . and iMaund et al.1 (1201 lh . The small red points 
represent the observed positions of Galactic RSG (iLevesque et al.ll2005l) . 



both the SN type and the estimated position in the HRD of the 
progenitor. 

One of the most important physical process governing the 
evolution of massive stars is mass loss. Combined with rota- 
tion, it changes the evolution of the star radically (see the re- 
view by Maeder & Mevnet 201 1). However, the mass-loss rates 
inferred observationally or theoretically remain uncertain. One 
of the most popular mass- l oss pr escriptions for the RSG phase 
is given by Ide Jaeer et al.l d 19881). and has been confir med ob- 
servationallv (ICrowfherl200UlMauron & Josselinll201 lb . On the 
other hand, some authors claim that the mass-l oss rates during 
the R SG phase cou ld have been underestim ated (Ivan Loon et al.1 
l200HlDavies et al.ll2008T: ISmith et al.ll2009l:lMoriva et al.ll201 ll) . 
Several authors have already studied the effects of such an in- 
creas ed mass-loss rate on the ev olution of non-rotating single 
stars. IVanbeveren et all (Il998albl) show that the RSG and WR 
lifetimes critically depe nd on the mass- l oss pr escription used 
during the RSG phase. lEldridge & Tout! (120041) studied partic- 
ularly the effects on the remnant mass and the final mass and 
ejected mass. They show that an increased mass-loss rate leads to 
lower remnant mass, and smaller hydrogen content in the ejecta. 



In lEkstrom et al.l (1201 lh . an increased mass-loss rate during 
the RSG phase was already used for stars with an initial mass 
M > 20 Mq during the RSG phase. In this letter, we extend this 
prescription to lower initial mass models, and explore system- 
atically the consequences of an increased mass-loss rate on the 
final position of the model in the HRD, as well as on the SN 
type produced at the end of the stellar life. In Sect. [2] we briefly 
recall the physics included in our models, particularly the mass- 
loss prescriptions used. In Sect. [3] we describe the effect of the 
increased mass-loss rate on the evolutionary tracks in the HRD. 
In Sect. 3] we discuss the consequences on the supernova pro- 
genitors and compare them with some recent observations. We 
present our conclusions in Sect. [5] 



2. Physics of the models 

The models computed for this work contain ex actly the same 
physic al ingredients as the models presented in Ekstrom et al.l 
(1201 lh . We mention here the main points of interest for this 
work: 

- The initial abundances are X — 0.720, Y = 0.266, and 
Z = 0.014, with a sola r mixture for the heavy elements 
dAsplund etal.ll2005L and lCunha etal]|2006l for the Ne abun- 
dance). This corresponds to the initial composition of the 
Sun. 

- The models are computed with the inclusion of the effects 
of the rotatio n, a ccording to the theoreti cal developments of 
IZahnl (Il992h and iMaeder & Zahnl (I 1998b . Vertical shear and 
meridional circulation are included in this approach. 

- In the mass domain studied in this paper , the mass-loss 
rates are taken from Ide Jager et al.l (JT988) for stars with 
log(T eff /K) > 3.7. For stars w ith log(T eff /K) ^ 3.7 , we 
use a linear fit from t he d ata of Sylvester et al.l d 19981) and 
Ivan Loon et all (1 19991) (see lCrowtherll200lh . 

- Due to the uncertainty on the mass-loss rate in the low-T e ff 
part of the HRD (see Sect. [TJ, we apply several prescrip- 
tions in this paper. We apply a multiplicative factor between 
three and ten to the standard mass-loss rate prescription de- 
scribed above, when the effective temperature of the star is 
lower than log(r e jf/K) < 3.8. The detailed list of the com- 
puted models is shown below (see Table[TJ. The exact phys- 
ical mechanism of this increased mass loss still has to be de- 
termined. One of the candidat es for such a mechanism is the 
pulsa tion-driven winds (Heger et al l 1 9971: lYoon & Cantiellol 
2010). Binary interactions c an also lead to an increa sed mass 
loss during the RSG phase (IVanbeveren et al.ll2007l) . 

To study the impact of an increased mass-loss rate during 
the RSG phase, we computed rotating stellar models of 12 and 
15 Mq, with an initial equatorial velocity given by v/v cl -i t = 0.4 
(which reproduces the observed mean velocity of B stars well 
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Table 1. List of the computed models. 



label 


Mini 


M 


label 


Mini 


M 




[Mq] 


[M -yr-'] 




[Mq] 


[Mq-vt- 1 ] 


12Std 


12 


Std 1 


15x3 


15 


Std-3 


12x10 


12 


Std • 10 


15x5 


15 


Std -5 


15Std 


15 


Std 1 


15x10 


15 


Std • 10 



Notes. (1) See text. 



during the main sequence (MS), see Ekstrom e t al.ll20i"il) . We 
performed computations with the standard and increased mass- 
loss rates. The list of the computed models is shown in Table [TJ 
Most of the models were computed up to the end of central car- 
bon burning. After that point, the central evolution of the star is 
supposed to be so quick that the surface properties (L, T^) are 
no longer modified. However, to verify this assessment, we com- 
puted two models up to the end of central neon burning (models 
15Std and 15x10). After this the star only has a few years before 
core collapse so the observable nature of the star will change 
only slightly. Luminosity and effective temperature of models 
computed up to the end of central Ne-burning only changed by 
less than 0.004 dex. Only one model, 15x5, was computed only 
up to the end of the central helium burning, and its final position 
should be considered a lower limit to the surface temperature 
and luminosity. 



3. Impact on the stellar evolution 

In the left-hand panel of Fig.Q] we show the evolutionary tracks 
of our models in the HRD. As a compa rison, we also plot th e 
track of the rotating 20 Mq model of Ekstrom et aT] (1201 ll) . 
which was also calculated with a mass loss increased by a factor 
of three during the RSG phase. The models computed with the 
standard mass-loss prescription (black and red tracks for the 12 
and 15 Mq models, respectively) follow a classical evolution- 
ary path, crossing the HRD after the MS, and then climbing 
the RSG branch at roughly constant T e e. The models with in- 
creased mass-loss rates (green, purple, and orange tracks for the 
mass-loss rates increased by factors of 3, 5, and 10, respectively) 
show the same behaviour. After they start to climb the super- 
giant branch, they all become hotter, going back in the bluer part 
of the HRD in a blue loop. They then come back in the yellow 
supergiant area, where they end their lives. We also see that the 
bluewards extension of the loop depends on the mass-loss rates: 
the higher the mass-loss rate, the hotter the star becomes during 
the loop (up to log(r eff /K) > 4). 

The mass-loss rates of the models are shown in Fig. [2] and 
the time evolution of their total mass in Fig. QO The models with 
a standard mass-loss prescription show a strong increase (1.5 
orders of magnitude) in the mass-loss rates as the star enters 
the RSG phase compared to the mass-loss rate during the MS. 
The models with increased mass loss reach strong mass-loss 
rates of around 10 [Mq yr~'] when they enter the RSG phase. 
During that phase, the mass decreases quickly. The models with 
increased mass loss then enter the blue loop, during which the 
mass-loss rates are lower. Because the duration of the loop is 
longer for stars with stronger mass-loss rates (up to ~ 50% of the 
total duration of the central helium burning for the models with 
mass-loss rates increased by a factor of 10), the final masses of 
the models with increased mass loss are comparable. 
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Fig. 2. Mass-loss rates of the models as a function of time, for 
the last million years of the stellar life. The labels near the curves 
indicate the corresponding model. 

4. Supernova progenitors 

4. 1 . Comparison with observed YSG supernova progenitors 

In the right-hand panel of Fig.Q] we show a zoom in the HRD, 
centred on the region where the end points of the evolutionary 
tracks lie. We also indicate the observed position of four re- 
cently discovered YSG SN progenitors: 1993J, 2008cn, 2009kr, 
and 201 ldh. All the suspected progenitors have a metallicity 
near solar, except S N 1993 J, whose m etallicity is estimated at 
roughly twice solar ( Sm artt et al.ll2002l) . Also these SNe were of 
the quite rare types IIL or lib (except SN 2008cn, classified as 
a type IIP SN). Classical evolutionary tracks of single star are 
unable to explain the position of these stars in the yellow area 
of the H RD. The binary channel is usually evoked to solve this 
probl e m dPodsiadlowski et al.1 19921 Il993blstancliffe & Eldrid"gel 
120091 IClaevs et al.l 1201 lh. H owever, in at least one case (SN 
201 ldhl. lMaund etal.l(l201lh indicate that no companion to the 
observed progenitor is detectable, making it more probable that 
this progenitor is a single stafl 

Looking at the end points of the tracks of the models with in- 
creased mass-loss rate, we see that they are in excellent positions 
to explain these YSG progenitors, covering the same T e s range 
as the observations. The progenitors of SN 201 ldh and 2009kr 
could thus be a 15 Mq star having encountered a stronger mass- 
loss rate during the RSG phase than currently adopted in stellar 
evolution codes. The slightly higher luminosity of the progen- 
itors of SN 1993 J and 2008cn indicates that they could have a 
slightly h igher initial mass t han our 15 Mq. As already men- 
tioned in Sma rtt et al.l (120091) . the final luminosity of our mod- 
els are roughly independent of the post-MS mass-loss history. 
Models with similar initial mass and rotation velocity end within 
a ±0.1 dex luminosity ranged The initial mass of an SN progen- 



2 This doe s not mean tha t this star was not a binary in the past. 
According to Eldrid ge~et al.1 d2011h . secondaries can be at the origin 
of up to 15% of all type lib, IIL and Iln SNe. 

3 However, varying the initial rot ation velocity also leads to a similar 
dispersion (see Ekstrom et al. 2 01 lh . 
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Fig. 3. Total mass of the models as a function of time. As a com- 
parison, the final mass of the CO-cores are 2.06 for model 12Std 
and 2.12 M for 12x10, and 2.87, 3.01, 2.67, and 2.74 Mq for 
the models 15Std, 15x3, 15x5, and 15x10, respectively. 




Fig. 4. Logarithm of the mass fraction of 'H, 4 He, 12 C, 14 N, and 
I6 at the surface as a function of the actual mass of the star, for 
the model 15Std (dashed line) and 15x10 (solid line). 



itor can thus be estimated using the luminosity of the progenitor 
alone. 

The surface abundances of the main chemical elements as a 
function of the actual mass of the star are shown in Fig. [4] For 
the model with standard mass-loss rates, we see that the surface 
chemical composition progressively evolves, revealing the prod- 
ucts of H-burning through the CNO-cycle: a decrease in the H, 
C, and O abundances and an increase in the He and N abun- 
dances. In the case of an increase in the RSG mass-loss rate by 
a factor of ten, the same trend is observed, but much stronger. 
In that case, the surface H abundance is strongly depleted, be- 
coming smaller than the He abundance. Measurements of the 
surface abundances indicating such trends would be a clue to 
strong mass loss during the RSG phase. 

As three of the four observed SNe discussed in this paper are 
of type IIL or lib, it is also important that our rotating single star 
models are compatible with such an explosion. As discussed in 
Heg er et al.l d2003). determining the supernova type on the basis 
of stellar models alone is not simple, as it depends strongly on 
the explosion properties. However, the absence of a plateau in 
the lightcurves of the types IIL and lib SNe indicates that the 
hydrogen envelope ejected during the e xplosion mu s t be s mall. 
We adopt here the sam e criterion as in iHeger et ail (120031) and 
lEldridge & Toutl (|2004|) to distinguish between type IIP and IIL 
SNe: it is a type IIL if the ejected mass of hydrogen is lower than 
2 Mq, otherwise it is a type IIP. 

The ejected hydrogen amounts during the explosion are 
3.8 Mq, 0.1 1 M Q , and 0.03 Mq for the models 15Std, 15x3, and 
15x10 respectively. The model with standard mass-loss prescrip- 
tion will thus produce a type IIP SN, and the other two a type IIL 
/ lib SN. An increased mass-loss rate during the RSG phase pro- 
duces rotating single-star models that fit both the position in the 
HRD and the SN type simultaneously for the cases we consider 
in this paper. Also the quite low hydrogen content of the models 
with increased mass loss at the end of their li ves could ex plain 
the lack of type IIP SNe from high-mass RSG (ISmartlll2009h . 



We can wonder to what extent an increased mass-loss rate 
is usual during the RSG phase. As already mentioned, there are 
theoretical and observational arguments for and against an in- 
creased mass-loss rate, making it very uncertain. In the right- 
hand panel of Fig. [T] the positions of Ga lactic RSG are indicated 
(small red circles. iLevesque et alj|2005l) . The time spent on the 
RSG branch depends on the mass-loss rate adopted. It corre- 
sponds to the whole central He-burning for the standard mass- 
loss prescription and to approximatively half of the same time 
for the models with mass loss increased by a factor of ten, be- 
fore the star evolves in the blue part of the HRD. Models with 
increased mass mass loss are thus able to reproduce the location 
of the RSG-br anch, as well as th e models with standard mass 
loss (see also|Ekstrom e t alJl20lH) . However, that the identified 
YSG SN progenitors exhibit a rather large dispersion in T e g, and 
on the other hand, there are a few type IIP SNe arising from stars 
with a n initial mass between ~ 15 - 18 Mq, which are probably 
RSG dSmarttl 12009). indicates that "standard" and "increased" 
mass-loss rates RSG coexist, leading to the observed variety of 
type II SNe and their progenitors. 

A way to determine the relative importance of each type of 
evolution would be to compare the relative number of type IIL/b 
SNe and type IIP SNe arising fro m stars with roughly the same 
initial mass. The actual data (e.g. Smartt et al. 2009) seem indi- 
cate that an increased mass-loss rate is relatively infrequent, due 
to the smaller number of type IIL/b SNe than type IIP. However, 
most of type IIP SNe arise from stars with initial mass lower than 
the ones considered in the present paper. If for some reason, an 
increased mass-loss rate is favoured for more massive stars (as 
expected, for example, in the case of the pulsation-driven winds 
discussed in lYoon & Cantiellol l2010h . the relative contribution 
of both mass-loss scenarios should depend on the initial mass, 
and the physical processes leading to the increased mass-loss 
rate could be more effective for stars of around 15 Mq than for 
smaller ones. 

One can also wonder whether the strong mass loss encoun- 
tered by our models leave observable imprints in the SN spec- 
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tra. The stellar winds during the RSG phase are slow winds, 
and might inte ract with the shock wave or the ejecta during the 
SN explosion. Filippenkg d 19971) mentions that type Iln SNe 
are SNe with a spectrum showing a strong interaction with a 
dense circumstellar medium, which indicates a strong mass-loss 
episode before the SN explosion. We cannot determine whether 
this should be the case here on the basis of our models. However, 
it would be interesting to find a observable signature of a strong 
mass loss in the SNe spectra. 

4.2. Low-luminosity Wolf-Rayet stars 

A particularly challenging question for the evolution of massive 
star is the existence of low-luminosity WC stars. They are not 
expla ined by current rotating single-star models dGeorgy et al.l 
1201 ll submitted), which produce much more luminous WC 
models. In this previous paper, we pointed out three solutions 
to this problem: 

- an increased mass loss during some part of the evolution of 
the most massive stars (during the MS and/or during the pre- 
vious WR phases), 

- an increased mass loss during the RSG phase of stars around 
15-20 Mq (the case explored in this work), 

- close binary evolution, where some mass transfer occurs be- 
tween the two components, allowing for the primary to lose 
its H-rich envelope. The effect on the evolution is similar to 
the previous item, instead that the mass loss is due to binary 
interactions (Eldrid ge et al.ll2008l) . 

The calculations performed in this work indicate that, even 
with a strong mass loss during the RSG phase, our models are 
still far from becoming WR stars. Even if the surface abundances 
show a depletion of hydrogen and an increase in helium, they do 
not fulfil the criteri on needed to be classified as a WR star (see 
iGeorgy et al J 12009). Moreover, it would be necessary to com- 
pletely remove the H-rich envelope of the star to become a WC 
star. In view of the structure of our models at the end of the stel- 
lar lifetime, it seems very unlikely that low-luminosity WC stars 
can arise from stars in the 15 - 20 Mq mass range through the 
single-star channel. 

5. Conclusions 

Motivated by observational and theoretical indications that the 
current prescriptions for the mass-loss rate of RSG could be un- 
derestimated and by the increasing number of YSG SN progen- 
itors, we performed a set of simulations of the evolution of ro- 
tating stars between 12 and 15 Mq, with various mass-loss rates 
during the RSG phase. 

We show that these kinds of models are able to reproduce 
the positions of the YSG progenitors in the HRD, fitting the ef- 
fective temperature well where these progenitors are observed, 
as well as the characteristics of the ejecta, leading to type IIL or 
lib SNe. Even if the exact physical mechanism leading to these 
high mass-loss rates still has to be determined, the consequent 
stellar evolution could thus be an alternative to the binary chan- 
nel, which is also able to reproduce the observed characteristics 
of the progenitor. 

Moreover, the low hydrogen content of the YSG progenitor 
will favour an explosion as a type IIL/IIb SN, giving a possi- 
ble explanation for the lack of massive RSG as progenitors of 
type IIP SNe. A more exhaustive study of the consequences of 
such an increased mass-loss rate will be addressed in a future 



paper, with a more extended set of mass, rotational velocities, 
and metallicities. 

Finally, the observed dispersion in effective temperature of 
the identified YSG progenitors, as well as the occurrence of mas- 
sive type IIP SNe, indicates that all the stars do not undergo the 
same mass-loss rate increase. However, a set of models with both 
various enhancements of the mass-loss rates and standard mass- 
loss rate is able to reproduce the whole range of observed pro- 
genitors. 
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